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We use the rigorous vector theory of weak photon localization for a semi-infinite medium 
composed of nonabsorbing Rayleigh scatterers to compute the full angular profile of the 
polarization opposition effect. The latter is caused by coherent backscattering of unpolarized 
incident light and accompanies the renowned backscattering intensity peak. 
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1. INTRODUCTION 


Coherent backscattering is a remarkable phenomenon caused by constructive interference of 
waves propagating along reciprocal multiple-scattering paths in a discrete random medium. 1-4 A 
well-known manifestation of coherent backscattering is an intensity peak centered at exactly the 
backscattering direction. The magnitude of the peak relative to the incoherent background 
depends on the polarization states of the incident and reflected beams and can be is large as a 
factor of 2 when the incident beam is fully circularly polarized, the intensity of the reflected 
beam is measured in the same-helicity channel, and the scattering particles are spherically 
symmetric. 3 

Another interesting and practically important manifestation of coherent backscattering can 
be observed when the incident beam is unpolarized. In this case the backscattering intensity 
peak is accompanied by a sharp asymmetric peak of negative polarization with a minimum 
centered at a very small phase angle (angle between the reflected beam and the backscattering 
direction). 6 This phenomenon was called the polarization opposition effect (POE) 6 and has the 
same physical origin as the so-called spatial anisotropy of the backscattering intensity peak in the 
case of fully linearly polarized incident beam. 7 

The calculation of the POE in Ref. 6 was based on the vector theory of coherent 
backscattering developed by Ozrin 8 for a semi-infinite medium composed of nonabsorbing 
Rayleigh scatterers (particle size is much smaller than the wavelength of the incident light). 
Although this theory is rigorous, the final solution was given in terms of asymptotic expressions 
valid in the limit of very small and very large phase angles. Therefore, the full angular profile of 
the POE, including the exact value and the angular position of the polarization minimum, 
remained unknown. 

In this paper, we revisit the problem of computing the full angular profile of the POE and, to 
this end, use the complete solution developed recently by Amic et al 9 We also compare our 
theoretical calculations with available laboratory data and briefly discuss remote-sensing 
ramifications of our results. 


2. THEORY 


The following analysis is largely based on formulas of Ref. 9. However, since we are interested 
in coherent backscattering of unpolarized incident light, it is more convenient to use the standard 
representation of the Stokes vector rather than the modified representation adopted in Ref. 9. 
This change necessitates a simple linear transformation of the 4x4 matrices given by Eqs. (2.42)- 
(2.45) of Ref. 9 as described in section 3 of Ref. 10. 

Consider a semi-infinite, homogeneous random medium composed of nonabsorbing 
Rayleigh scatterers. We will use the standard set of Stokes parameters /, Q, U, and V to define 
the polarization state of the incident and reflected beams relative to their respective meridional 
planes (the planes through the beams and the normal to the boundary of the medium). 1011 
Assume that the medium is illuminated by a parallel beam of light incident perpendicularly to the 
boundary and characterized by Stokes parameters /o, Oo, Uq, and Vq, where /o is the incident 
energy flux per unit area perpendicular to the flux. The Stokes parameters ot the beam reflected 
at a phase angle ot in the meridional plane of the incident beam are given by 
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where S(a) is the Stokes reflection matrix. 10 The reflection matrix can be decomposed as 1 ' 


S(a) = S'(a) + S £ (a) + S c (a) , ( 2 ) 

where S'(a) is the first-order-scattering component, S A (a) is the contribution ot all ladder 
diagrams of scattering orders 2 and higher, and S ( (ct) is the contribution of all cyclical diagrams 
of scattering orders 2 and higher. Restricting the analysis to small ct, introducing the so-called 
angular parameter q = kla , where k is the wavenumber and / is the photon transport mean tree 



path in the scattering medium, and performing the previously mentioned matrix transformation, 
we derive from formulas of Ref. 9 
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Since the angular functions y ^(q) can be computed numerically as described in detail by Amic 

et al „ 9 Eqs. (2H 10) can be used to calculate the full Stokes reflection matrix for any q. If, 
furthermore, the transport mean free path / is known, the reflection matrix can be expressed as a 
function of the phase angle a rather than the dimensionless angular parameter q. Note that the 
matrices S^O) and S c (0) given by Eqs. (4)— (10) fully comply with the general relationships 
derived for arbitrary macroscopically isotropic and symmetric scattering media and given by 
Eqs. (18) and (20)-(23) of Ref. 5. 

Assuming unpolarized incident light (Qo — Uo~ Vq — 0), we define the intensity enhancement 
factor as the ratio of the total reflected intensity to the incoherent background value, 


,(_) S|',(0) + S|y0) + S|S(<?) (11 ) 

’’ SiVo+s/.iO) 

Equations (3)— <5) show that only the first two Stokes parameters of the reflected light can be 
nonzero. Therefore, the degree of linear polarization can be defined as minus the ratio of the 
total reflected second Stokes parameter to the total reflected intensity: 

P( a) = -Q^l = giztoj . (12) 

1(a) ShW + SkW + Sfa) 

Figure 1 shows both quantities as functions of the dimensionless angular parameter q. The upper 
panel demonstrates the renowned intensity peak centered at exactly the opposition. The 
amplitude of the peak is ^(0) * 1.5368 and its half-width at half-maximum is qr * 0.597. The 
former number agrees with the previous calculation reported in Table 2 of Ref. 5. whereas the 
latter number should replace the value 0.56 reported previously in Ref. 6 and derived using the 
asymptotic formulas of Ref. 8. Thus the relationship between the half-width at half-maximum of 
the backscattering intensity peak and the photon transport mean free path for the case of 
conservative Rayleigh scattering and unpolarized normal illumination is given by 
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The lower panel of Fig. 1 shows that the reflected polarization is zero at opposition, which is 
a consequence of azimuthal symmetry in the case of normal illumination and unpolarized 
incident light. However, with increasing q . polarization becomes negative, rapidly grows in 
absolute value, and reaches its minimal value P mm ~ -2.765% at a reflection direction very close 
to opposition ( qp » 1 .68). The peak of negative polarization is highly asymmetric so that the 
half-minimal polarization value -1.383% is first reached at q Py \ « 0.498, which is even smaller 
than the value q^ « 0.597 corresponding to the half-width at half-maximum of the backscattering 
intensity peak, and then at a much larger qp y2 ~ 7.10. This unusual behavior of polarization at 
near backscattering angles was called in Ref. 6 the polarization opposition effect. To the best of 
our knowledge, the numbers q ; * 0.597, P,™,, =» -2.765%, qp « 1.68, q Py \ * 0.498 and q Py2 « 7.10 
provide the first exact metrics of the photometric and polarization opposition effects caused by 
coherent backscattering of unpolarized incident light by a half-space of nonabsorbing Rayleigh 
scatterers. 

3. DISCUSSION 

Because lasers are the most frequently used sources of illumination and usually generate linearly 
or circularly polarized light, explicit laboratory demonstrations of the POE have been extremely 
rare. Apparently the first laboratory observation of the POE was made by Lyot as long ago as in 
the 1920s, 12 although the physical explanation of this phenomenon was, of course, unknown at 
that time. Figure 2 shows Lyot’s polarization measurements for a particulate surface obtained by 
burning a tape of magnesium under a glass plate until the deposit on the plate was completely 
opaque. Lyot described the observed phase curve of polarization as “puzzling” and attributed it 
to the very small size of magnesia grains. Unfortunately, he did not measure the actual size of 
the grains and their packing density and, thus, did not provide the information necessary to 
compute the photon mean free path /. Furthermore, the minimal measured polarization value is 
only -1.1 1%, compared to the theoretical value P m in ~ -2.765% computed for Rayleigh particles. 
However, assuming that the latter difference is explained by the finite particle size in the Lyot’s 
experiment and multiplying the theoretical polarization by a factor of 0.4, and also assuming that 
the actual kl was close to a realistic value of 132, we were able to almost perfectly reproduce the 
angular profile of the mesaured polarization up to phase angles of about 15° (solid curve in Fig. 
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2). At larger phase angles the assumption of phase-angle-independent single-scattering and 
ladder contributions to the total reflection matrix [Eqs. (3) and (4)] is no longer valid and causes 
a significant deviation of the theoretical curve from the measurements. 

The POE may have also been observed in more recent laboratory measurements by Geake 
and Geake 13 for fine alumina powders, although the smallest phase angle in their experiments 
was 2° and may be too large to make definitive conclusions. These measurements may suggest 
that the POE weakens and ultimately disappears with increasing particle size, which would be 
consistent with the similar disappearance of the spatial anisotropy of the backscattering intensity 
peak in the case of linearly polarized incident beam observed by van Albada et al: 

Sharp peaks of negative polarization at phase angles well less than a degree were observed 
for Saturn’s rings and the Galilean satellites of Jupiter/ ' 14 Because these peaks were distinctly 
different from the symmetric, nearly parabolic negative polarization branch commonly observed 
for many solar system bodies, 1 ' they w r ere interpreted in terms of the POE produced by coherent 
backscattering of sunlight by a surface layer composed primarily of microscopic ice particles." 14 
The fact that the same objects exhibit comparably narrow intensity peaks centered at exactly the 
opposition 16 provides a strong support for this interpretation. 

As was emphasized previously, our present analysis is based on rigorous formulas of Ref. 9 
and provides, for the first time, the exact metrics of the POE in the case of conservative Rayleigh 
scattering. Since the vector theory of coherent backscattering is extremely complicated and 
contains only a few rigorous results, 5 ’ 8 ’ 9 we believe that our calculations can be useful in 
checking the accuracy of approximate approaches that are often based on multiple untested 
assumptions and have not been validated versus exact theories. 17 Furthermore, it appears that 
laboratory and remote-sensing measurements of polarization at nearly zero phase angles can be 
an informative characterization and diagnostic tool in various physical and geophysical 
applications. 
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FIGURE CAPTIONS 


Fig. 1. Theoretical angular profiles of the intensity enhancement factor and the degree of linear 
polarization of the reflected light for a half-space of nonabsorbing Rayleigh particles illuminated 
by an unpolarized beam of light incident normally to the boundary of the scattering medium. 

Fig. 2. Polarization measurements for a particulate surface composed of microscopic magnesia 
particles (dots) and best-fit theoretical computations (solid curve). 
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